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1.  Introduction 

A  wind  atlas  is  an  important  tool  for  optimal  exploitation  of 
wind  resources  [1,2].  The  offshore  wind  atlas  of  the  Canary  Islands 
[3,4]  is  a  geographic  chart  that  gives  information  about  three 
important  wind  features:  mean  wind  speed,  wind  direction  and 
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wind  potential  at  different  locations  around  the  islands.  In  the 
following  sections  an  explanation  has  been  included  of  the  steps 
that  have  been  followed  to  compile  data  and  establish  the  charts 
for  this  atlas,  and  graphs  are  given  to  reveal  the  results  obtained. 

2.  Stages  of  the  work 

The  atlas  has  been  developed  in  the  following  six  stages: 

1.  Definition  of  the  research  area. 
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2.  Data  collection  from  different  anemometer  stations. 

3.  Obtaining  of  maps  from  numerical  models  and  satellites. 

4.  Statistical  data  processing. 

5.  Calculation  and  elaboration  of  the  offshore  wind  atlas. 

6.  Data  extrapolation. 

In  the  following  sections  the  six  previous  steps  are  explained  in 
detail. 

3.  Definition  of  the  research  area 

The  first  step  for  the  elaboration  of  a  wind  atlas  (regardless  of 
whether  it  is  offshore,  inshore  or  nearshore)  consists  of  defining 
the  working  area  limits,  which  constitute  what  has  been  denoted 
as  working  domain  in  this  work.  In  this  case  the  limits  of  the 
geographic  area  are  given  by  the  parallels  of  northern  latitude 
29  °  55'00"  and  26°  41 '00",  and  the  meridians  of  western 
longitude  18  °  50'00"  and  12  °  38'00". 

This  means  the  studied  area  comes  to  212,760  km2,  and 
contains  the  7  main  islands,  i.e.,  El  Hierro,  La  Palma,  La  Gomera, 
Tenerife,  Gran  Canaria,  Luerteventura  and  Lanzarote,  and  the  6 
smaller  islands,  i.e.,  Alegranza,  La  Graciosa,  Montana  Clara, 
Lobos,  Roque  del  Este,  Roque  del  Oeste,  together  with  a  small 
part  of  the  African  coast,  belonging  to  Western  Sahara 
[5,6].  Geographic  information  handled  in  this  work  has  been 
given  by  SITCAN,  which  is  the  Sistema  de  Information  Territorial 
de  Canarias  (Territorial  Information  System  of  the  Canary 
Islands). 

4.  Data  collection  from  different  anemometer  stations 

Anemometer  stations  for  data  collection  must  be  sought 
according  to  the  following  criteria: 

1.  These  stations  have  to  measure  the  sought  after  data,  of  course, 
such  as  wind  speed  and  direction  in  this  case. 

2.  It  is  important  to  have  data  for  3  years,  at  least. 

3.  They  should  be  located  as  close  as  possible  to  the  area  being 
studied. 

4.  They  have  to  be  rigorously  controlled  when  acquiring,  storing, 
correcting  and  verifying  data. 

Data  from  40  weather  stations  have  been  handled  for  the 
elaboration  of  this  atlas.  Their  locations  on  the  different  islands  are 
the  following:  2  stations  are  on  El  Hierro,  2  on  La  Palma,  2  on  La 
Gomera,  15  on  Tenerife,  11  on  Gran  Canaria,  3  on  Luerteventura 
and  3  on  Lanzarote.  2  stations  are  buoys  floating  in  the  Atlantic 
Ocean,  one  of  them  to  the  south  of  the  island  of  Tenerife  and  the 
other  one  at  the  northwest  of  Gran  Canaria.  So,  the  quantity  of 
acquired  data  comes  to  2,102,400,  i.e.,  2  data  (wind  speed  and 
direction)  per  hour  during  3  years  (24  x  365  hours  per  year)  for 
the  40  stations. 

5.  Obtaining  of  maps  from  numerical  models  and  satellites 

Comparative  studies  have  shown  that  correct  simulation  of  a 
phenomenon  depends  on  the  specific  configuration  and  options 
used  for  its  modelling.  It  is  not  so  important  to  work  with  a  given 
numerical  model  for  weather  prediction,  but  rather  to  choose  an 
adequate  parametrization  and  configuration  [7].  A  bibliography 
review  has  been  made  into  different  modelling  configurations, 
interaction  techniques,  resolution,  and  model  application  [8-11]. 
Configurations  of  operational  models  in  several  meteorological 
institutions  have  been  studied,  and  also  some  research  works 
where  meteorological  models  have  been  applied  to  the  analysis  of 
weather  phenomena. 


The  final  conclusion  has  been  that  a  good  solution  is  to  apply  a 
technique  known  as  ‘one  way  nesting’,  that  is  to  say,  interaction  in 
one  direction  [12].  This  technique  groups  together  results  obtained 
from: 

1.  Satellite  mapping  (QuickSCAT). 

2.  Limited  area  modelling  (HIRLAM). 

3.  Mesoscalar  modelling  (MM5). 

Another  requirement  for  this  work  is  to  have  very  powerful 
hardware  for  calculations.  This  can  be  substituted  by  requesting 
data  from  those  international  institutions  that  are  able  to  provide 
them. 

As  for  satellite  mapping  we  should  comment  that  data  have 
been  provided  by  the  NASA  (US  National  Aeronautics  and  Space 
Administration)  QuickSCAT  satellite  [13,14],  which  contains  a 
microwave  radar  sensor  that  can  measure  reflexion  or  dispersion 
produced  as  it  explores  the  Earth’s  surface.  This  radar  was 
specifically  developed  for  measuring  wind  speed  and  direction 
near  the  ocean  surface  and  has  three  subsystems:  SES  (electronic 
subsystem),  SAS  (antenna  subsystem),  and  CDES  (command  and 
data  subsystem). 

The  SES  is  the  core  of  the  module  and  contains  a  transmitter,  a 
receiver  and  a  processor  for  numerical  signals.  It  generates  and 
sends  radiofrequency  waves  to  the  SAS,  which  sends  the  signal  to 
the  Earth  by  means  of  pulses.  These  pulses  reach  the  ocean  surface 
and  experience  an  echo  effect,  known  as  backscatter.  This 
backscatter  is  redirected  by  the  SAS  to  the  SES,  and  this  system 
converts  the  received  signal  to  facilitate  its  processing.  The  CDES  is 
basically  a  computer  with  the  software  needed  to  operate  the 
whole  system  correctly.  One  of  the  operations  the  CDES  has  to 
carry  out  is  to  calculate,  with  minimum  error,  the  position  of  the 
signal  on  the  Earth’s  surface.  The  CDES  also  stores  the  temperature 
of  the  system  and  operating  currents  and  voltages  to  supervise  the 
state  of  the  device. 

Data  given  by  this  satellite  are  particularly  important  for 
drawing  up  the  offshore  wind  atlas  as  they  allow  us  to  know  not 
only  the  interaction  between  earth  and  ocean  but  also  about 
ocean  streams,  together  with  their  effects  on  global  climate 
behavior. 

In  Fig.  1  an  example  is  shown  of  a  map  of  ocean  winds  at  a 
height  of  10  m.  Data  have  been  processed  by  NOAA/NESDIS 
(National  Oceanic  and  Atmospheric  Administration/National 
Environmental  Satellite  Data  and  Information  Service)  [15],  from 
data  obtained  by  NASA/JPL  (Jet  Propulsion  Laboratory),  by  means  of 
QuicSCAT.  These  images  have  a  resolution  of  20  km  and  have  been 
adapted  to  a  grid  made  specifically  for  the  Canary  Islands  offshore 
wind  atlas.  This  grid  contains  31  columns  and  21  rows,  which 
results  in  651  smaller  rectangles,  not  all  of  them  completed. 

The  fact  that  images  are  captured  once  per  day  should  been 
taken  into  account.  Data  corresponding  to  4  years  have  been 
requested  (period  2003-2006),  which  means  a  total  of  1,900,920 
data.  The  mean  value  for  every  rectangle  is  calculated  from  all 
these  data. 

By  applying  a  ‘one  way  nesting’  technique,  these  data  will  serve 
as  a  reference  for  correcting  the  HIRLAM  model  [16,17].  The  goal  of 
the  correction  is  to  avoid  excessive  deviations  in  the  final  results, 
and  consists  of  establishing  maximum  and  minimum  values  at 
115%  and  85%  of  the  calculated  mean  value.  So,  values  outside 
these  limits  are  corrected  and  restricted  to  them. 

The  HIRLAM  model,  in  operation  in  the  Spanish  Weather 
Agency  AEMET  (formerly  known  as  INM),  has  the  following  main 
features: 

1.  Two  possible  horizontal  resolutions:  0.16  °  and  0.05  °. 

2.  Vertical  resolution:  40  hybrid  levels. 
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Fig.  1.  Wind  speed  and  wind  direction  obtained  from  the  QuickSCAT  satellite. 


3.  Mesh  resolution:  9  km. 

4.  First  guess:  predictions  based  on  registers  taken  6  h  before 
with  the  same  model. 

5.  Analysis: 

(a)  Height: 

i.  Data  assimilation:  intermittent  every  6  h  with  a  data  window 
of  -3  to  +2  h. 

ii.  Variables:  temperature,  wind,  relative  humidity. 

iii.  Method:  three-dimensional  variational. 

(b)  Surface: 

i.  Variables:  ocean  water  temperature,  pressure  at  ocean  level, 
wind  speed  at  10  m  altitude  and  temperature  at  2  m. 

ii.  Method:  univariant. 

6.  Initialization:  by  digital  filtering. 

7.  Prediction: 

(a)  Dependent  variables:  pressure  at  surface  level,  temperature, 
wind  and  specific  humidity. 

(b)  Vertical  representation:  finite  differences,  energy  conserva¬ 
tion,  angular  moment. 

(c)  Horizontal  representation:  finite  differences,  Arakawa  C 
grid. 

(d)  Integration  time:  Leapfrog,  semi-implicit  with  time  step 
equal  to  2  or  4  min. 

(e)  Horizontal  diffusion:  implicit  4th  order. 

(f)  Vertical  diffusion:  limit  layer. 

(g)  Dynamic:  Semilagrangian. 

An  image  extracted  from  the  map  generated  by  the  AEMET  is 
shown  in  Fig.  2. 

From  maps  received  from  AEMET  [18],  a  mesh  has  been  built 
with  a  resolution  of  9  km,  consisting  of  66  rectangles  in  west-east 
direction  (horizontal),  and  40  rectangles  in  north-south  direction 
(vertical),  which  gives  2640  computable  rectangles. 


It  has  been  applied  in  the  running  of  the  HIRLAM  model  for  the 
5-year  period  2002-2006.  This  model  presents  data  for  4  daily 
passings  (at00:00, 06:00, 12:00  and  18:00).  By  calculating  the  total 
number  of  data  on  the  basis  of  the  product  explained  earlier,  the 
number  of  total  data  handled  is  38,544,000,  where  wind  speed  and 
direction  are  included.  The  mean  values  are  obtained  to  calculate 
one  rectangle. 

Similar  to  what  happens  with  satellite  images,  the  technique  of 
interaction  in  one  direction  is  employed,  such  that  the  data  of  the 
9  km  domain  will  serve  as  reference  for  the  correction  of  the  MM5 
model  [19].  Such  a  correction  consists  of  restricting  values  above 
115%  or  below  85%  of  the  value  calculated  by  HIRLAM,  as 
mentioned  before. 

Focusing  on  the  MM5  model,  whose  data  are  obtained  by  the 
ITER  ( Instituto  Tecnologico  de  Energias  Removables )  [20]  and  the  RIU 


Wind  Tested  to  10  m  (knots) 
10  07  2006  at  00  hours  UTC 


Fig.  2.  Image  by  HIRLAM  model  extracted  from  AEMET. 
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( Red  Iberica  de  Usuarios  de  MM5)  [21  ],  the  most  relevant  aspects  of 

the  version  used  (3.6)  are  the  following: 

•  Capacity  for  multiple  nesting  in  both  directions  between 
domains,  which  facilitates  the  analysis  of  atmospheric  phenom¬ 
ena  for  different  time  scales,  and  the  design  of  very  high 
resolution  predictions. 

•  Formulation  of  a  non-hydrostatic  dynamic,  which  allows  the 
effective  use  of  the  model  to  represent  phenomena  of  very  small 
size  in  km. 

•  Data  adaptation  for  using  on  multiple  platforms  and  multitask¬ 
ing  execution  in  parallel  computing. 

•  Automatic  initialization  with  different  sources  of  meteorological 
analysis  and  observations,  including  the  ability  to  assimilate 
four-dimensional  data. 

•  Introduction  of  the  most  modern  and  realistic  parametrization 
schemes  for  physical  processes  with  regard  to  atmospheric 
radiation,  micro-physics  of  clouds  and  precipitations,  cumulus 
convection,  turbulence  and  energy  flows  and  moment  over  the 
Earth’s  surface. 


speed  in  m/s  and  wind  direction  at  a  given  altitude  Zg.  The 
reference  for  wind  direction  is  the  east  direction.  North  direction 
means  90  °  and  so  on. 

So,  data  from  anemometer  stations  and  those  from  prediction 
models  are  given  in  the  same  reference  frame  and  measurement 
units,  in  the  International  System. 

7.  Calculation  and  elaboration  of  the  offshore  wind  atlas 

The  offshore  wind  atlas  for  the  Canary  Islands  has  been  defined 
on  the  basis  of  mean  wind  speed  and  mean  direction  at  10  m 
height.  The  values  employed  are  the  result  of  the  interaction 
between  three  domains  of  the  three  prediction  models  (satellite 
mapping,  HIRLAM  model  and  MM5  one).  The  results  are  given  for  a 
resolution  of  3  km.  These  values  are  corrected  with  data  from  the 
anemometer  stations.  The  coastal  data  for  the  islands  are  the  result 
of  an  interpolation  between  the  values  of  prediction  models  with  a 
weight  of  65%,  and  the  value  taken  at  the  nearest  anemometer, 
with  a  weight  of  35%.  At  6  km  of  the  coast  this  correction  is 
negligible. 


The  data  for  initializing  the  model  come  from  the  global  model 
GFS  (Global  Forecast  System)  [22,23]  for  the  initial  contour 
conditions.  A  relatively  high  horizontal  resolution  was  chosen  for 
MM5  to  be  able  to  show  meteorological  phenomena  that  affect  the 
islands  with  enough  detail. 

As  a  test  of  the  results  given  by  MM5  model,  an  assessment  has 
been  made  for  some  days  with  a  modern  model,  WRF  (Weather 
Research  and  Forecasting  model)  [24],  and  the  conclusion  was  that 
the  results  were  almost  identical. 

The  numerical  values  of  wind  speed  and  direction  given  by 
MM5  belong  to  the  period  2004-2006,  with  a  resolution  of  3  km, 
which  involves  a  mesh  of  1 97  rectangles  in  west-east  direction  and 
120  in  north-south  direction.  This  means  23,640  rectangles. 

The  MM5  model  uses  data  from  8  daily  passings  (at  00:00, 
03:00,  06:00,  09:00,  12:00,  15:00,  18:00  and  21:00).  Multiplying 
the  number  of  rectangles  by  the  number  of  years,  the  days  per  year, 
the  number  of  passings  and  2  (because  there  are  two  variables, 
speed  and  direction),  a  total  amount  of  414,172,800  data  are 
obtained.  The  mean  value  is  what  is  taken  into  account. 

As  a  summary,  the  total  number  of  data  used  with  the  three 
models  is  454,617,720.  Due  to  this,  a  data  sheet  has  been  used  for 
each  model  and  year,  using  Open  Office  Calc.  The  results  of  the 
maps  are  refreshed  automatically  when  there  are  new  data. 

6.  Statistical  data  processing 

Due  to  the  large  amount  of  data  to  be  processed,  that  is  to  say, 
454,617,720  from  the  prediction  models,  and  2,102,400  from  the 
anemometer  measurements,  a  very  exhaustive  control  of  data 
collection  is  required. 

6.2.  Prediction  model  data  analysis 

Once  these  data  have  been  introduced  in  the  different  data 
sheets,  an  algorithm  has  to  be  run  to  search  for  values  greater  than 
35  m/s  or  lower  than  0  m/s,  and  erroneous  values.  In  this  last  case 
they  are  substituted  by  assumed  valid  values  by  means  of  an 
interpolation  of  the  nearest  values  in  the  data  register. 

Data  from  the  different  prediction  models  are  compared  with 
those  of  the  verifying  equation  in  order  to  check  if  they  are  identical. 


7.2.  Mean  wind  speed 

Mean  wind  speeds  at  a  height  of  10  m  can  be  seen  in  Fig.  3.  This 
map  shows  interesting  information  about  the  best  locations 
according  to  wind  speed.  The  areas  with  best  features  seem  to  be 
the  south  and  southeast  part  of  the  island  of  Gran  Canaria,  the  area 
between  Tenerife  and  west-northwest  of  Gran  Canaria,  southeast 
of  the  island  of  El  Hierro,  northeast  of  La  Palma,  and  northwest  of  La 
Gomera.  In  all  these  locations  mean  wind  speeds  can  reach  values 
in  the  interval  of  8-10  m/s. 

In  the  case  of  Lanzarote  and  Fuerteventura  the  mean  wind 
speeds  are  in  the  interval  of  4-6  m/s,  and  the  worst  areas  from  the 
point  of  view  of  wind  speed  are  located  in  the  southwest  of 
Tenerife  and  La  Palma,  with  values  between  3  and  4  m/s. 

7.2.  Wind  direction 

On  the  other  hand  the  mean  wind  direction  maps  can  be  seen  in 
Figs.  4-6. 

According  to  the  8  main  directions  and  to  the  reference 
explained  (0  °  for  east  direction),  the  following  notation  has  been 
included  in  the  maps: 

•  East:  from  337.5  °  to  22.5  °,  symbol-*. 

•  Northeast:  from  22.5  °  to  67.5  °,  symbol 

•  North:  from  67.5  °  to  112.5  °,  symbol  |. 

•  Northwest:  from  112.5  °  to  157.5  °,  symbol  \. 

•  West:  from  157.5  °  to  202.5  °,  symbol  . 

•  Southwest:  from  202.5  °  to  247.5  °,  symbol 

•  South:  from  247.5  °  to  292.5  °,  symbol  j. 

•  Southeast:  from  292.5  °  to  337.5  °,  symbol  \. 

By  looking  at  the  maps  we  can  say  that  the  wind  comes  mainly 
from  the  north,  although  there  are  certain  areas  where  other 
directions  dominate,  such  as  northwest,  northeast  or  east. 

Although  the  mean  value  has  been  represented  on  the  maps,  an 
interesting  subject  is  to  check  if  this  value  is  actually  representa¬ 
tive,  for  which  two  statistical  values  have  been  calculated: 

The  first  one  is  a  measurement  of  dispersion,  and  is  given  by  the 
standard  deviation,  defined  in  (1): 


6.2.  Anemometer  station  measurements  test 

Wind  speed  data  have  been  taken  from  anemometer  stations 
existing  inside  the  working  domain.  Each  station  provides  the  wind 


E  (*/-*) 2 


n 


(1) 
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Fig.  3.  Mean  wind  speed. 


where: 

•  Xi  is  a  given  speed  or  direction  value. 

•  n  is  the  number  of  speed  or  direction  samples. 

•  X  is  the  mean  speed  or  mean  direction  value,  which  is  calculated 

as  X  =  J2i=\  xi/n- 

The  mean  wind  speed  seems  to  be  a  more  reliable  value  than 
the  direction,  because  its  standard  deviation  always  remains  at 


values  in  the  interval  of  3-4  m/s,  whereas  the  standard  deviation 
for  direction  reaches  values  of  100  °. 

There  is  also  a  measurement  of  shape  and  concentration,  given 
by  the  Fisher  asymmetry  coefficient  (2): 

X>-*)3-/,- 

? -  (2) 


■*■  v  •*  -ir  « 


rvr'rwvr'T'ri 


r  rr  v  r  i 


1.  V  te  i. 


M  1»  A  A  A 


*  •*  A  A  A 


A**  ******  AAAA* 
X  A ’A  A 


A  A  A  A  A  AAA-AAAAAA 

aaaaaaaaaaaaaa 

:;:a:aa;aa;a:. 

1  AAilAUAAAU  . 
4’4Ml'T4''fT'r,r’r'rvrv 

m  a  i  i  n  a  u  i  a  : 
AAAAAAAAAAAAAA 


AAAAAAAAAAAAAA* 


A4AAAAAAAAAA 

AAAAAAAAAAAA 

AAAAAAAAAAA. 

AAAAAAAi-AAAA 


AAA AAAAAAAAAAA 
A AAAAAAAAAAA A A 
AAAAAA-A  +  A-A-A-AA^ 
AAAAAAAAAAAAAA 


•  *«AAA  AAAAAAAAi  A*. 
AAAAAAAAAAAAAA4AAA 
...  A  A  A  A  A  A  A  A  i  A  A  A  ..  . 
A  ;  «  A  A  AAA  AAAAAAA*»  A ■ 
AAAAAAAAAAAA**; 
AAAA AAAAAAAAAAAAAA 
AAAAAAAAAAAAAAAAAA 


AAAAAAAAAAA. 

AAAAAAAAAAAA 

AAAAAAAAAAA* 

AAAAAAAAAAAA 

AAAAAAAAAAA. 


AAAAAAAAAAAA. 
AA'AA' A  AAAAAAAA' 
AAAAAAAAAAAAa*. 
AAAAAAAAAAAA; 


AAAAAAAAAAA 
AAAAAAAAAAAAAA 
AAAAAAAAAAAA.. 
AAAAAAAAAAA**. 
AAAAAAAAAAAAAA 
AAAAAAAAAAAAAA 
AAAAAAAAAAAAAA 
AAAAAAAAAAAAAA 
AAAAAAAAAAAAAA 
AAAAAAAAAAA.. . 
AAAAAAAAAAA*;* 
AAAAAAAAAAAAAA 
AAAAAAAAAAA. * • 
AAA AAAAAAAAAA* 
A^AAA;AiAAA;;; 
AAAAAAAAAAA.;. 
AAAAAAAAAAAAAA 
AAAAAAAAAAAAAA 

aaaiaaaaaaa;;. 

AAAAAAAAAAAAAA 

♦  A  A  A  A AAAA 
AAA AAAAAA 
AAAAAAAAA 
A  A  AAA AAAA 

AAAAAAAAA 


AAAAAAAAAAA 

AAAA*AAAAAA 

AAAAAAAAAAA 

AAAAAAAAAAA 

AAAAAAAAAAA 


A  AAAA 


.  . 


.  . 


.  . 


;  ;  ;  •  •  «  v  ;  •  m  ;  ;  ;  ;  'w 

*  *  *w  ^  v.  *w  *;  *»**•*»  '» 

'<  *.  A  A  A 

; *  «  *  *  *  A  *  AAAAAA 
A^\<^*^AAAAAAAA 

•  .  '.  A  A  A  1  .  A  A  A  A  A  A 


AAAAAAAAAA 
AAAAAAAAAA 
AAAAAAAAAA 
A  A  ;  A  AAAAAA 
AAAAAAAAAA 


AAAAAAAAAA 

AAAAAAAAAA 

AAAAAAAAAA 

AAAAAAAAAA 

AAAAAAAAAA 

AAAAAAAAAA 


r AAA A A A AAA 
►AAAAAAAAA 
AAAAAAAAA 
AAAAAAAAA 
.AAAAAAAAA 


AAA 

AAA 


A  A 


A  A 


AAA 


A  A 


^  ^  ^  ^  ^  ^  ^  ^  ^ 
AAAAAAAAAAAAA 
A  A  A  ^  ^  ^  ^  ^  ^  ^ 

AAAAAAAAAAAAA 

^  ^  ^  ^  ^  ^  ^ 


AAAA 


A  VL 

^  ^ 

AAAAAAAA 

>A'  v  4*  v*  *A*  ^ 

AAAAAAAA 

^  ^  ^  vjy  ^  ^ 

^  ^  ^  ^ 

^  ^  ^  ^  ^  ^  ^  ^ 


AAAA 

AAAAAAAAAAAAAAAAAAAAAAAA 

■  ^y  ^y  ^y  \_^y  ^y  «^y  -^r 

A.  A  A.  A  A  ^  ^  ^  ^  A  AAAAAAAAAAAA 

•X^ir  AAAAAAAAAAAA 


^  ^  ^  ^  ^  ^  ^  ^  ^ 


^  ^  ^  ^  ^  ^ 

X  X  X*  X  X  >1*  X  X  X  X  X 

XXXXXXXX 

'X'*  ^y»  ^y>  X*  X** 

\  X  X  'A'  X  A  X  >A 

X  x  ^  ^  ^  >^y  ^  ^ 


\.|y  ^y  v|y  ^ 

X  X  X  X  X  X  X  X  X  X  X  X  X 

^  .^y  ^i‘ 


X  X  X  X  X  * 


Fig.  4.  Mean  wind  direction  -  western  area. 


Fig.  5.  Mean  wind  direction  -  middle  area. 
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Fig.  6.  Mean  wind  direction  -  eastern  area. 


where: 

•  Xj,  n,  X  and  S  mean  the  same  as  in  (1). 

•  fj  is  the  frequency  of  each  value. 


According  to  this  coefficient,  the  distribution  can  be  considered 
symmetrical  if  y=  0.  It  can  be  said  that  both  mean  speeds  and 
directions  tend  to  a  symmetrical  distribution,  although  in  the  case 
of  mean  speeds,  values  close  to  2  or  -2  are  also  reached. 

A  new  and  different  set  of  three  maps  have  been  drawn  up  and 
represented  in  Figs.  7-9,  where  the  dominant  direction  of  the  wind 
has  been  represented. 

The  conclusion  is  that  the  northeast  direction  is  dominant  for 
85%  of  the  domains,  although  north  and  northwest  are  also 
dominant  directions  in  some  cases. 

The  assumption  has  been  made  that  there  is  no  need  to  estimate 
mean  directions  and  mean  speeds  over  a  long  term  period  by 
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Fig.  8.  Dominant  wind  direction  -  middle  area. 


means  of  methods  such  as  Corotis  [25]  or  Box-Jenkins  [26].  This  is 
because  the  fact  has  been  assumed  that,  in  the  least  favorable 
situation,  a  series  of  3  years  measurements  is  available  for  study. 
According  to  research  works  by  Corotis,  this  means  that  there  is  a 
90%  confidence  level  for  the  mean  value  of  direction  or  speed  over 
the  next  20  years. 

As  for  number  of  operating  hours,  the  Weibull  distribution  has 
been  chosen  as  representative  of  this  kind  of  parameter.  This  is  a 
well  known  distribution,  given  by  two  coefficients,  scale  (C),  and 
shape  (/<),  that  can  be  put  in  relationship  with  the  mean  value  and 
standard  deviation  by  means  of  the  moments  method.  The  Weibull 
distribution  can  be  expressed  using  its  probability  distribution 
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Fig.  7.  Dominant  wind  direction  -  western  area. 


Fig.  9.  Dominant  wind  direction  -  eastern  area. 
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function  according  to  (3). 

P{Xi  <x0)  =  1  -  e~(Xo/C)*  (3) 

where  P(x*  <  x0)  denotes  the  probability  that  x*  <  x0. 

The  values  of  I< and  C  can  be  calculated  according  to  (4)  and  (5): 

/  c\  — 1.086 

K  -  (I)  (4> 


c  =  x 


Tl+K- 


where  T  is  the  Gamma  function  of  Legendre. 

When  trying  to  check  the  fitting  level  between  the  approxi¬ 
mated  function  and  the  observed  distribution,  the  following 
aspects  can  cause  differences: 


1.  Fluctuations  during  the  sampling  process. 

2.  The  observed  distribution  does  not  correspond  to  the  law  as  was 
assumed. 


In  order  to  establish  the  fitting  level,  a  Smirnov-Kolmogorov 
test  [27]  has  been  carried  out,  by  defining  a  hypothesis 
H0 :  F(x)  =  F0(x),  to  compare  both  distributions.  H0  is  rejected 
when  there  is  a  significant  difference  between  both  distributions. 
In  mathematical  terms,  the  test  is  defined  by  means  of  (6)  and  (7): 

Dn  =  max  |Sn(x)  -  F0(x)|  (6) 

[  0  x  <  xk 

Sn (x)  =  <  k/n  xk  <x<xM  (7) 

[  1  X  >  Xn 

The  distribution  of  Dn  is  independent  of  the  proposed  model 
under  H0,  with  the  result  that  it  can  be  assessed  exclusively  as  a 
function  of  the  sample  size.  The  error,  a,  can  be  established  as  a 
function  of  a  critical  region  given  by  (8): 

p(Dn>7n)=a  (8) 

Kolmogorov-Smirnov  test  is  not  used  in  all  those  cases  where 
the  observations  are  inherently  quantitative  due  to  the  ambiguities 
that  can  appear  when  putting  the  observations  in  order. 

Once  the  method  and  test  that  have  been  used  for  the 
calculation  of  the  operating  hour  maps  have  been  commented 
on,  it  becomes  necessary  to  explain  that  two  kinds  of  estimators 
have  been  defined:  the  first  one  by  stretches  (around,  under  and 
above  the  mean  value);  and  a  different,  full  type  one.  Several  maps 
can  been  elaborated  for  these  options,  although  they  are  not  all 
included  in  this  paper. 

7.2.1.  Estimation  of  operating  hours  around  the  mean  value  of  the 
calculated  speed 

The  number  of  hours  during  which  a  given  generator  will  be 
working  in  an  interval  close  to  the  mean  speed  is  calculated  on  the 
mean  value  map.  Such  an  interval  is  determined  by  upper  and 
lower  limits  of  5%  around  the  mean  value. 

Losses  are  estimated  as  2%  (due  to  non-operating  hours  for 
maintenance  and  the  like),  which,  over  1  year,  results  in  175.2  h. 

7.2.2.  Estimation  of  operating  hours  under  the  calculated  mean  speed 
The  number  of  hours  during  which  the  generator  will  be 

operating  under  the  calculated  mean  speed  is  estimated  according 
to  the  mean  speed  map.  The  interval  that  is  computed  for  the 
estimation  begins  in  the  cut-in  value  (that  we  are  considering 
equal  to  3  m/s),  and  goes  up  to  the  95%  of  the  mean  speed  at  that 
point. 


It  is  also  assumed  that  global  losses  reach  values  of  4%,  i.e., 
350.4  h. 

7.2.3.  Estimation  of  operating  hours  above  the  calculated  mean  speed 
The  number  of  hours  during  which  a  given  generator  will 

operate  above  the  mean  speed  is  estimated.  The  interval  is 
computed  from  105%  of  the  mean  speed  and  30  m/s. 

A  value  of  4%  is  also  considered  for  the  losses. 

By  analyzing  the  previous  maps,  the  conclusion  can  be  that  the 
total  operating  hours  are  in  the  interval  of  1500-2000. 

7.2.4.  Estimation  of  operating  hours  grouped  by  the  wind  potential 
index 

The  operating  hours  valid  for  wind  generation  at  nominal 
power  are  estimated  taking  into  account  the  design  at  each 
domain.  As  the  nominal  power  is  very  variable  according  to  the 
manufacturer  and  features  of  each  machine,  an  interval  given  by 
the  most  usual  speeds  has  been  used.  When  dealing  with  offshore 
wind  turbines,  these  are  7,  8,  9,  10,  11  and  12  m/s. 

The  wind  potential  index  (WPI)  allows  us  to  classify  the 
different  areas  according  to  the  number  of  equivalent  operating 
hours  (eh)  per  year  valid  for  wind  generation. 

Four  categories  have  been  established: 

•  WPI=  1  — >  e/i  <  2900  h 

•  WPI  =  2  2900  h  <  eh  <  3300  h 

•  WPI  =  3  ->  3300  h  <  eh  <  3700  h 

•  WPl  =  4  ->eh>  3700  h 


According  to  the  previous  map,  it  can  be  observed  that  almost 
90%  of  the  calculated  domain  results  in  WPI  =  3.  There  are  three 
areas  with  WPI  =  4,  one  of  them  located  to  the  northwest  of  La 
Gomera,  a  second  one  to  the  west  of  Gran  Canaria  and  south  of 
Tenerife  and  a  third  one  to  the  east  of  Gran  Canaria.  The  area  with 
lowest  value  of  WPI  is  to  the  southwest  part  of  Tenerife,  including  a 
part  of  the  coast  of  La  Gomera. 

As  a  summary,  comparing  data  from  the  mean  wind  speed  map 
and  data  from  the  operating  hour  estimation  maps,  the  conclusion 
can  be  reached  that  the  highest  wind  speeds  (interval  from  8  to 
10  m/s)  and  the  greatest  number  of  operating  hours  can  be  reached 
in  three  areas:  The  eastern  part  of  Gran  Canaria,  with  an 
approximate  surface  of  135  km2,  the  southeastern  part  of  Tenerife 
and  the  western  part  of  Gran  Canaria,  with  a  surface  of 
approximately  200  km2,  and  the  northeastern  part  of  La  Gomera, 
with  an  equivalent  surface  of  100  km2. 

8.  Data  extrapolation 

One  of  the  most  meaningful  phenomena  when  exploiting  wind 
energy  is  wind  shear.  Due  to  the  friction  of  the  air  mass  with  the 
Earth  surface,  wind  speed  decreases  from  an  undisturbed  value  at 
great  height  to  a  value  of  0  at  ground  level.  The  Earth  limit  layer 
consists  of  a  number  of  layers,  each  of  them  governed  by  a  different 
set  of  flux  parameters.  However,  the  layer  closest  to  the  surface  and 
the  Ekmann  layers  are  the  most  interesting  for  structure  design 

[28] .  The  surface  layer,  from  the  ground  to  approximately  100  m 
height,  is  the  region  where  the  variation  of  the  shear  tension  can  be 
neglected,  and  this  is  the  area  where  the  wind  turbines  are  located. 

Hiester  and  Pennell  have  pointed  out  that  only  measurements 
taken  at  the  height  of  the  hub  of  a  wind  turbine  will  provide 
enough  accuracy  for  a  reliable  calculation  of  the  wind  resource 

[29] .  However,  for  preliminary  estimations,  it  is  possible  to  reduce 
the  cost  of  taking  measurements  at  great  heights  by  using 
measurements  taken  at  a  reference  height,  and  then  extrapolating 
them. 
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Generally,  a  neutral  atmosphere  is  considered,  and  a  value  of 
roughness,  z0,  is  estimated.  Then,  (9)  is  employed  to  estimate  the 
mean  wind  speed  at  a  height  z,  knowing  the  wind  speed  at  a  height 
h. 


V(z) 


v(h)  ■ 


(9)  is  a  well  known  equation,  frequently  employed  in  energy 
analysis  in  Europe.  For  use  in  the  offshore  wind  atlas  of  the  Canary 
Islands,  data  are  extrapolated  to  4  typical  hub  heights  (40,  60,  80 
and  100  m).  The  results  can  be  represented  on  several  maps. 

By  means  of  an  exhaustive  analysis  of  all  the  maps  obtained 
after  the  extrapolation  at  these  different  heights,  we  can  accept 
that  there  is  no  great  variation  between  the  results  calculated 
initially  for  10  m  and  those  from  this  extrapolation,  due  to  the  low 
level  of  roughness  of  the  ocean. 


9.  Estimating  offshore  wind  power  generating  capacity  in  the 
Canary  Islands 


As  a  general  balance  we  have  analyzed  the  offshore  wind  power 
that  can  be  obtained  in  the  Islands  with  wind  turbine  generator 
technology  existing  at  present. 

On  the  basis  of  the  maps  obtained  for  the  mean  value  and  the 
estimation  of  operating  hours,  the  studied  domain  is  classified 
according  to  the  wind  potential,  from  which  the  wind  power 
capability  in  offshore  locations  has  been  calculated.  The  following 
hypotheses  have  also  been  considered: 


•  A  very  general  estimate  has  been  made  for  the  marine  areas 
located  inside  Spanish  territorial  waters. 

•  For  this  purpose,  wind  turbine  generators  between  1 .5  and  5  MW 
have  been  employed. 

•  The  wind  turbine  generator  distribution  on  the  terrain  has  been 
assumed  on  the  basis  of  the  dimensions  of  the  machines 


Table  2 

Offshore  wind  energy  that  can  be  obtained  with  different  types  of  wind  turbine 
generators. 


Island 

Wind  energy  (GWh)  I 

Wind  energy  (GWh)  II 

La  Palma 

608.27 

2742.57 

El  Hierro 

380.62 

1216.54 

La  Gomera 

380.56 

1905.55 

Tenerife 

1171.80 

3871.33 

Gran  Canaria 

1511.33 

6873.83 

Fuerteventura 

1572.78 

4817.50 

Lanzarote 

1170.50 

3298.86 

Total 

6795.86 

24526.18 

It  is  also  very  important  to  understand  that  Table  1  shows  wind 
power  potential  in  offshore  areas,  but  that  we  do  not  consider  here 
electrical  network  capability  for  integrating  such  an  amount  of 
power,  a  task  which  is  beyond  the  remit  of  this  paper.  A  deeper 
analysis  of  such  data  allows  us  to  think  about  the  fact  that, 
although  the  islands  of  Fuerteventura  and  Lanzarote  present 
higher  values  of  possible  power  to  be  installed,  this  does  not  mean 
that  wind  speeds  are  higher  there. 

Using  previously  described  criteria  we  can  venture  to  show  what 
production  by  means  of  offshore  wind  parks  would  be.  Assuming  the 
described  Weibull  distribution  for  this  task,  in  Table  2  we  can  see  the 
results  obtained,  also  for  both  kinds  of  wind  generator. 

A  comparison  between  the  24,526  GWh  of  estimated  offshore 
wind  production  and  the  trend  in  forecasted  demand  for  the 
period  2008-2015  for  the  Islands,  the  conclusion  is  that  offshore 
wind  production  would  be  able  to  cover  the  demand  not  only  in 
2015,  but  perhaps  also  until  2030.  This  would  not  be  a  real 
situation  due  to  power  system  penetration  limits  and  other 
factors. 

10.  Conclusions 


employed  during  the  study,  and  also  the  distribution  of  wind 
directions  has  been  taken  into  account,  for  each  area,  according 
to  data  obtained  from  the  Dominant  Wind  Direction  Map  given 
in  Figs.  7-9. 

•  The  results  of  power  capacity  are  given  according  to  the  kind  of 
foundation  the  generator  has,  because  this  has  influence  on  the 
depth. 

•  Locations  with  WPI  <  3  are  not  accepted.  In  all  these  marine 
domains  there  are  areas  that  could  be  environmentally 
sensitivity  and  have  been  declared  protected,  these  include 
natural  areas,  bird  conservation  areas,  marine  life  areas,  and  so 
on.  These  areas  have  to  be  excluded  from  wind  park  installation. 

In  Table  1  results  have  been  given  of  the  total  offshore  wind 
power  that  could  be  installed  for  two  kinds  of  wind  turbine 
generators,  i.e.,  type  I  (with  a  foundation  fixed  at  the  bottom  of  the 
sea),  and  type  II  (floating). 


Table  1 

Offshore  wind  power  that  can  be  installed  with  different  types  of  wind  turbine 
generators. 


Island 

Wind  power  (MW)  I 

Wind  power  (MW)  II 

La  Palma 

171 

771 

El  Hierro 

107 

342 

La  Gomera 

139 

696 

Tenerife 

428 

1414 

Gran  Canaria 

332 

1510 

Fuerteventura 

846 

2592 

Lanzarote 

589 

1660 

Total 

2612 

8985 

In  this  paper,  an  offshore  wind  atlas  for  the  Canary  Islands  is 
presented.  Such  an  atlas  is  based  on  three  combined  models,  by 
computing  a  wide  data  set.  Three  maps  have  been  drawn  up,  one  for 
Mean  Wind  Speed,  one  for  Mean  Wind  Direction  and  one  for 
Dominant  Wind  Direction.  Studying  them  allows  us  to  conclude  that 
a  great  offshore  wind  power  resource  is  available  in  the  islands.  The 
estimate  of  offshore  wind  power  generating  capability  also  indicates 
a  high  value. 
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